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Abstract 

We report a search for the decay ip(2S) — > 7?7 C (2S') in a sample of 25.9 x 10 6 ip(2S) events 
collected with the CLEO-c detector. No signals are observed in any of the 11 exclusive r] c (2S) 
decay modes studied, or in their sum. Product branching fraction upper limits are determined as 
a function of F[r] c (2S)} for the 11 individual modes. 
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I. INTRODUCTION 

The first radially excited S-wave spin singlet state in the charmonium system, t] c (2S), 
was observed by the Belle Collaboration in the decay process B ± — >■ K ± r] c (2S),T] c (2S) — > 
KgK ± ir Zf It was confirmed by the CLEO [2[ and BaBar [3[ Collaborations in the two- 
photon fusion process e + e" — > e + e~(77),77 — > rj c (2S) — > K^K^n^ and by the BaBar 
Collaboration in the double-charmonium production process e + e~ — > J/i/jcc [4j. These 
observations, which give an average mass M[r] c (2S)} = 3638 ± 4 MeV/c 2 , are inconsistent 
with a previous measurement of M[r] c (2S)] = 3594 ± 5 MeV/c 2 |5] based on an inclusive 
measurement of the decay process ip(2S) — > , yr] c (2S). By measuring the inclusive photon 
spectrum in a sample of 1.6 x 10 6 ip{2S) decays collected with the CLEO III detector, the 
CLEO Collaboration set an upper limit of B(ip(2S) ->• ir] c {2S)) < 0.2% at 90% confidence 
level (C.L.) for an Vc {2S) mass of 3594 MeV/c 2 @. 

Theoretical predictions for the branching fraction of 'M2S) — > , ~fri c (2S) based on potential 
model calculations fall in a range of (0.1 - 6.2) x 10~ 4 [7] for M[r] c (2S)} = 3638 MeV/c 2 . A 
phenomenological prediction, based on assuming that the matrix element governing ip{2S) —¥ 
7?7 C (25') is the same as that for J/ip — >■ •yrj^lS), is given by 

BmS) -> 7^(25)) = ^ B(J/^ -> 1V c(lS)), (1) 

fcj/^ t V(25) 

where k^^s) is the photon energy for the ip(2S) — >■ •yr] c (2S) [J /if) — >■ 7^(15')] transition, 

IV(2S) [r w ] is the V(2S) [J/V>] Ml width, and #(J/V> -»■ JVcO-S)) = (1.72 ± 0.25)% is the 
weighted average of the value listed by the Particle Data Group (PDG) jsj and a recent 
CLEO measurement Using the PDG values for k^s), kj/^, T^s), and T j/^, leads to a 
prediction of B(i)(2S) -)• 7^(25)) = (3.9 ± 1.1) x 10~ 4 . 

In this paper we describe a search for rj c (2S) production through ip{2S) —> r yrj c {2S) 
using a sample of 25.9 x 10 6 ip{2S) decays collected with the CLEO-c detector. We at- 
tempt to fully reconstruct the r] c (2S) in 11 exclusive decay modes: KKtt (composed of the 
decay modes K^K^ir* and K + K~tt°), 2(7r+7r-), 3(7r+7r-), K+K'tt+tt-, K+R-n+n-n , 
K + K~2(tt + tt~), KgK 7t^tt + tt~ , TT + TT~f], K + K~f], it + tt~t]', and 7i + 7i^r] c ( IS) , where the 
r) c (lS) is reconstructed in KKtt, 2(tt + tt~), and K + K~tt + tt~ decays. For a specific de- 
cay mode T) C (2S) — > X, the yield of events, N sig , in a sample of ij){2S) decays, N^s)i is 
given by 

N sig ^eN il{2S) B(^(2S)^ 1 r ]c (2S))B( Vc (2S)^X), (2) 

where e is the efficiency for fully reconstructing the r] c (2S) hadronic decay and the ~50 
MeV transition photon. The main experimental challenge is background from low-energy 
photons. By searching for exclusive decays to specific final states, it may be possible to 
observe r) c (2S) decays through modes other than KKtt, determine the product branching 
fractions in Eq. and obtain improved measurements of rj c (2S) properties. 



II. DATA SAMPLE AND EVENT SELECTION 

We use an e + e~ annihilation data sample with an integrated luminosity of 51.8 pb _1 taken 
at the ip(2S) mass, \/s = 3.686 GeV. The data were produced with the symmetric electron- 
position beams delivered by the Cornell Electron Storage Ring (CESR) and collected by the 



CLEO-c detector |10J. CLEO-c is an approximately cylindrically symmetric detector that 
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provides a solid angle coverage of 93%. The charged particle tracking system, consisting of 
a six-layer wire vertex detector (ZD) and a 47-layer wire drift chamber (DR), provides a 
momentum resolution of 0.6% for tracks with transverse momenta of 1 GeV/c. An electro- 
magnetic calorimeter (CC) consisting of 7784 cesium iodide crystals detects electromagnetic 
showers with an energy resolution for photons of 2.2% at E 7 = 1 GeV and ~ 5% at 100 
MeV. Charged particle identification (PID) information is obtained by measuring ionization 
energy loss (dE/dx) in the DR and with a Ring Imaging Cherenkov (RICH) detector. The 
RICH detector is located between the DR and CC and covers | cos6*| < 0.83, where 9 is 
defined with respect to the positron beam. All of the detector components described here 
reside within a 1.0 T magnetic field aligned with the beam axis. 

Reconstructed events are required to have the appropriate number of charged tracks for 
the exclusive process being investigated and therefore to have zero net charge. Charged 
tracks not associated with a Kg decay are required to have an impact parameter within 
5 mm of the e + e~ annihilation interaction point (IP) and within 5 cm of the IP along the 
beam axis. The charged pions used to form Kg candidates are constrained to a common 
vertex, which is required to be displaced from the IP by at least 3 standard deviations as 
determined from the net momentum of the 7r + 7r~ pair, and are required to have an invariant 
mass within 10 MeV/c 2 of the K s mass. 

PID information from dE/dx and the RICH detector is combined to discriminate be- 
tween pions and kaons. Separation is achieved by a requirement on the variable Axn = 
X 2 dE/dx( n ) ~ X 2 dE /dx( K ) - 2 In Crich{^) + 2 In C RICH (K), where XdE/dx(i) is the number of 
standard deviations of separation between the measured dE/dx and the mean expectation 
for a particular particle hypothesis, and Crich^) is the likelihood for a particular particle 
hypothesis using information from the RICH detector. We require kaon candidates to have 
Akit > 0, while pion candidates are required to have Ak-k < 0. If there is no information 
for the track from the RICH detector, the number of Cherenkov photons associated with 
the track is less than 3, or the momentum of the track is less than 700 MeV/c, then only 
dE/dx information is used. If a decay mode includes two charged kaons, then only one kaon 
is required to pass the A Kn criterion. 

Transition photon candidates are required to be detected in the region of the CC that gives 
the best performance (|cos#| < 0.81, "barrel"), to have a lateral shower shape consistent 
with that of a photon, to be not associated with a charged track traversing the CC, and 
to have a minimum energy of 30 MeV. The rj candidates are reconstructed in the rj — > 77 
and rj — > 7r + 7r~7T° decay modes, with the 7r + 7r _ 7r° invariant mass required to be within 10 
MeV/c 2 of the nominal rj mass [M(r/)]. The rj' is reconstructed from the decay process 
rf — > 7r + 7r~?7 with 77 — > 77, with the 7r + 7r — 77 invariant mass required to be within 10 MeV/c 2 
of M(rj'). Photon pairs forming a 7r° (rj) candidate are selected from both the barrel and 
endcap (0.85 < | cos#| < 0.93) regions of the CC. They are required to have a two-photon 
invariant mass within 3 standard deviations of the nominal mass, approximately ±18 (±36) 
MeV/c 2 for the ir° (rj), and are kinematically constrained to the 7r° (rj) mass for subsequent 
event reconstruction. 

Monte Carlo (MC) simulations, which have been extensively tested with independent 
data samples, are used to determine detector efficiencies and to study backgrounds. The 
MC samples are generated with EVTGEN [llj and a GEANT-based [l2[ detector simulation. 
Radiation emitted from charged particles, i.e., final state radiation (FSR), is simulated with 
the PHOTOS package [l3| • For signal MC samples, the generated angular distribution 
of the vector to vector-pseudoscalar ip(2S) — > ^r] c (2S) transition is 1 + cos 2 9, while the 
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T] C (2S) is decayed according to phase space. MC samples consisting of 259 x 10 6 generic 
ip(2S) decays (10 times the data size) and a "continuum" sample of e + e~ — >■ qq (q = u, d, s) 
events, consisting of an integrated luminosity of 259 pb _1 (5 times the data size), are used to 
study possible backgrounds. The generic ip{2S) MC sample is gene rated using the available 
branching fractions for the ij){2S\ y r j, J/ip, and r] c (lS) decays . with unmeasured decay 
modes simulated by JETSET (l5jj . The continuum sample is generated using JETSET 
models and has been validated with a data sample consisting of 20.6 pb _1 collected at 
yfs = 3.671 GeV. 

Transitions from the if)(2S) resonance to other low- lying charmonium states are a poten- 
tially large background for ip(2S) — > •yr] c (2S). In particular, decays to the J/ip via 7r + 7r~ 
and 7] transitions and transition photons from ip(2S) — > 'JXcj need to be suppressed. Selec- 
tion criteria to suppress these decays were designed and efficiencies and background- rejection 
fractions determined with signal and background MC samples, respectively. In what follows, 
efficiency loss and background rejection are determined from the differences in the yields 
with and without the requirement being investigated. 

For the hadronic final states K + K^-k + -k" , K + K~7t + 7t~tt , KgK ± n T 7T + 7i^ , 3(7r + 7r~), and 
K + K~2(7r + 7T~), events are rejected if any 7r + 7r~ pair originating from the IP has a recoil 
mass within 20 MeV/c 2 of M(J/ip) or if the invariant mass of the other hadrons is within 
30 MeV/c 2 of M(J/ip). Efficiency losses for K + K~tt + tt~ are 0.1% for both criteria, and 
backgrounds are reduced by 10% and 1% for the 7r + 7r~ recoil and invariant mass criteria, 
respectively. The efficiency loss for K + K~7T + 7T~7r is 0.4% (0.6%) for the 7r + 7r~ recoil (in- 
variant) mass criterion, while the background is reduced by 26% (4%). Efficiency losses for 
KgK ± n T n + 7Y^ are 1.0% for both criteria, and backgrounds are reduced by 57% and 34% for 
the 7r + 7r~ recoil and invariant mass criteria, respectively. Efficiency losses for 3(7r + 7r~) and 
K + K~2(ii + 'k~) are 5.5% and backgrounds are reduced by two-thirds for the 7r + 7r~ recoil 
mass criterion while, for the invariant mass criterion, the efficiency loss is 5.6% [6.3%] and 
the background is reduced by 33% [51%] for Z(ix + tj-) [K + K~ 2{tt + tt-)]. 

In order to suppress the much more abundant i[ ! (2S) — > 7r + 7r~ J/ip, J/ip — > £ + £~(£ = e, fi) 
decays in the 2(7r + 7r~) final state, events are rejected if the recoil mass of any 7t + tt~ pair 
originating from the IP has a value greater than M(J/ip) — 30 MeV/c 2 = 3067 MeV/c 2 . The 
efficiency loss for this criterion is 0.7%, while it reduces the background by 96%. 

To suppress ip{2S) — > T]J/ip decays, events are rejected if the 7] recoil mass is within 
40 MeV/c 2 of M(J/ip) for tt+tt^, 7] 77 and within 20 MeV/c 2 of M(J/ip) for n+n-r], 
f] — > 7r + 7r~7r°. Efficiency losses are 0.5% for both rj decays, while backgrounds are reduced 
by two-thirds. 

For the K + K~it° final state, events are rejected if the recoil mass determined from the 
higher energy photon used in forming the 7r° candidate is within 20 MeV/c 2 of M{\c2) or 
M(xd), or within 30 MeV/c 2 of M(x c0 ). The efficiency loss is 4%, while the background is 
reduced by 38%. The same photon recoil mass requirements are used for the lower energy 
photon in rj — > 77 decays for the 7r + 7r — 77 and K + K~rj final states. The efficiency loss is 24% 
(27%) for 7r + 7r~?7 (K + K~rj), while the background is reduced by 83% (71%). 

The invariant mass of the hadronic decay, M inv , for all modes is required to be between 
M[ip(2S)) and 100 MeV below it, i.e., AM = M[ip(2S))-M inv with < AM < 100 MeV/c 2 . 
Requiring AM > MeV/c 2 rejects events with a direct ip{2S) decay combined with a low- 
energy shower, while requiring AM < 100 MeV/c 2 rejects hadronic decays of the Xc2 state. 
The efficiency loss is largest for modes with only two charged tracks (6.6% for K + K~tt°, 
2.7% for K + K~i],i] — > 77, and 2.4% for 7r + 7r — 77, 77 — > 77), while it is less than 1.7% for all 
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other modes. The background rejection ranges from 21% for the K + K~2(tt + 7t~) mode to 
67% for the 2(7r + 7r~) mode. 

Kinematic fitting is used to optimize signal detection and reject background. The sum of 
the four-momenta of the reconstructed hadronic decay and the transition photon candidate 
is constrained to the initial ip{2S) four-momentum. The requirement on the \ 2 P er degree 
of freedom (x 2 /d.o.f.) for this total event fit is optimized mode by mode by evaluating 
the figure of merit S 2 /(S + B). The accepted signal (S) is determined by processing a 
sample of signal MC events that was generated with an assumed branching fraction of 
B{ip(2S) — > 7^(25*)) = 2.6 x 10 -4 [l6[ (smaller than our current phenomenological estimate) 
and the arbitrary assumption that the branching fraction for each rj c (2S) decay to light 
hadrons is 1%. The generic ip{2S) and continuum background MC samples, scaled to our 
data sample size, are used to compute the corresponding background (B). The xVd.o.f. 
requirements derived from this study are listed in Table HI In addition, the reconstructed 
particles originating from the IP are constrained to a common vertex, and the y 2 per degree 
of freedom of this vertex fit is required to be less than 10 for all modes. 

TABLE I: Mode-dependent full event fit x 2 /d.o.f. selection criteria and overall signal efficiencies 
(e) for the r] c (2S) decaying into light hadrons. Efficiencies include statistical uncertainties and 
constituent decay mode branching fractions 0]. We have assumed T[r] c (2S)] = 14 MeV/c 2 . 



Channel 




X 2 /d.o.f. 


e(%) 






<3.5 


14.09 ±0.10 


K+K-ir° 




<4.0 


17.55 ±0.14 


KKir 






7.63 ± 0.04 


2(vr+vr-) 




<4.5 


20.48 ±0.16 


3(vr+vr-) 




<5.0 


14.22 ±0.14 


K+K--K+-K- 




<4.0 


19.50 ±0.15 






<2.5 


8.68 ±0.11 


K+K-2(tt + 


TP") 


<4.0 


9.93 ±0.11 






<4.0 


7.84 ± 0.09 


ir + ir~ri, r] — > 


- 77 


<2.0 


4.03 ± 0.04 


ir + ir~ri, r] — > 


- 7T + 7r~7r° 


<3.0 


1.65 ±0.02 








5.68 ± 0.05 


K + K~r/, r] - 


-> 77 


<3.5 


4.55 ± 0.05 


K + K~r/, r] - 




<5.0 


1.92 ±0.02 


K+K~r] 






6.48 ± 0.05 






<3.0 


1.42 ± 0.02 



Additional selection criteria have been developed for suppression of low-energy shower 
backgrounds. These showers are associated with bremsstrahlung radiation emitted from 
charged pions in the reconstructed hadronic decays (FSR) and showers created from nu- 
clear reactions of charged pions and kaons in the CC ("split-off" showers). The x 2 /d.o.f. 
requirement for the total event fit suppresses some split-off showers, but it does not provide 
effective suppression of FSR since the energy momentum is balanced in a fully reconstructed 
hadronic decay with FSR. FSR can be suppressed by requiring that the opening angle be- 
tween a charged pion at the IP and the transition photon candidate be greater than some 
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value. Split-off showers can be suppressed by requiring the transition photon candidate to 
be some distance away from a charged track entering the CC. The specific selection criteria 
are optimized using the same S 2 / (S + B) procedure described above. 

While all decay modes are evaluated for additional background shower suppression, FSR 
suppression is found to be useful only for the 2(7r + 7r~) and 7r + 7r — 77, 77 — > 7r + 7r~7r° hadronic 
final states, for which the angle between a charged pion and candidate photon is required to 
be greater than 0.376 rad. The efficiency loss is 13.0% (11.1%) for the 2(7r + 7r~) (n + n^ri, 77 — > 
7r + 7r~7T°) mode, while the background is reduced by 41% (48%). Split-off suppression is 
only applied to the K + K^-k + -k^ and K + K~7r° hadronic final states, for which the distance 
between the charged track and candidate photon shower is required to be greater than 45 
and 35 cm, respectively. The efficiency loss is 7.7% (2.2%) for the K + K~7t + tt~ (K + K~ir° ) 
mode, while the background is reduced by 32% (18%). 



III. YIELD DETERMINATION PROCEDURE AND CROSS-CHECKS 

Searches for the ip(2S) —> •yr] c (2S) transition are performed by studying the measured 
shower energy of the transition photon candidate without adjustment from the total event 
kinematic fit. Signal yields are determined by performing a binned log-likelihood fit of the CC 
shower energy distribution with a Breit-Wigner function convoluted with a MC-determined 
detector resolution function for the signal shape and a background shape composed of the 
events from the background MC samples that pass the event selection criteria. The mean 
and width of the Breit-Wigner function are fixed to E 1 = 48 MeV and T 7 = 14 MeV j8j. 
Any monochromatic shower energy distribution reconstructed in the CC has a low-side tail 
caused by losses sustained in interactions prior to entering the CC and from leakage outside 



the Csl crystals. For that reason, the Crystal Ball function [17J is used to parameterize the 
detector resolution, with parameters determined from the signal MC samples. 

The procedure for determining the ij){2S) — > 7^(2^) yields has been studied and tested 
with two closely related processes. Reconstructed i/j(2S) — > r yXc2,Xc2 X decays are 
used to test the signal fitting procedure. To assess the reliability of the MC samples for 
determining the shape of the background shower energy distribution, we investigate the 
process ip{2S) — > 7T + 7T~J/tp with the J/tp decaying to 2(tt + tt~), K + K^-k + -k^ , K + K~tt°, 
or KgK ± 7f T . These modes have mixtures of final state hadrons very similar to our signal 
channels and no additional photons, so the calorimeter response should closely resemble the 
backgrounds in the rj c (2S) signal region. More details on these studies are provided in Ref. 
13. 



For the J/ip decay study, the previously described event selection criteria for the 3(7r + 7r~), 
K + 'K~2(7t + 7t~), K + K~7T + 7T~7r°, and KgK ± n T 7i + 7i^ hadronic final states are applied with 
the exceptions that the 7r + 7r~ recoil mass and J/ip hadronic decay suppression criteria are not 
applied and that at least one tt + tt~ pair is required to have a recoil mass within 20 MeV/c 2 
of M(J/ip). Three different background shapes were studied: a first-order polynomial, the 
energy distribution of showers from the background MC samples that pass the tt + it~J/iI) 
selection criteria (one free parameter for the normalization), and events from the same 
background MC samples partitioned into separate distributions for showers identified as 
being split-off showers and for all others (two free parameters, the normalization of each 
distribution). Figure [TJ shows the fits of the measured shower energy distributions with 
the backgrounds predicted by the MC with one free parameter and no special treatment of 
split-off showers. The xVd.o.f. are 25.1/20, 11.7/20, 11.7/20, and 39.7/20 for the J/rp ->• 
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FIG. 1: Distributions of measured shower energy for the decay modes ip(2S) —> tt + tt~J/iP, (a) 
J/ip -> 2(7r + vr-), (b) J/i/j -> K+K-7T+TT-, (c) J/ip -> K+K-7T , and (d)J/ip -> i^l?^. 
The points are data and the solid histograms are the background MC distributions with floating 
normalization. 



2(7r + 7r~), K + K~7t + tt~ , K + K~7c°, and KgK ± ir T decays, respectively. The background MC 
samples are found to adequately reproduce the behavior observed in data, with no clear 
improvement when split-off showers are treated separately. Therefore, we use the shower 
energy distributions from the background MC samples with a single normalization parameter 
in fitting the r] c (2S) signal region. 

For the ip(2S) —> 7Xc2, Xc2 — > X study, the ip{2S) — > r yr) c (2S),r] c (2S) — > X event selection 
criteria are applied, except that the transition photon is required to be in the energy range 
of 90 — 145 MeV and no AM requirement is applied. Figure [2] shows the shower energy 
distributions for the 2(7r + 7r~) and K + K~tt + 7t~7t° final states. The transition photon signal 
is well fitted, validating the MC-determined resolution function. Table [III lists the number 
of signal events observed and, for comparison, the number of events expected for our ip(2S) 



sample using branching fractions from the PDG [14|. The efficiencies are determined from 



signal MC samples for each hadronic final state, where the generated angular distribution of 
the ip(2S) —> 7Xc2 decay is 1 + ^ cos 2 9 (which assumes a pure El transition jl9j), while the 
Xc2 is decayed according to phase space. The energy resolution for the transition photon, 
determined from the signal MC samples, is ~6.2 MeV. Since the substructure of the Xc2 
decays and the systematic uncertainties of these measurements are not evaluated, these 
yields are presented only as a cross-check of the yield determination procedure and not as 
measurements of the Xc2 decays. 

The numbers of observed events listed in Table [H are consistent with the PDG lli 
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FIG. 2: Distributions of measured shower energy for the decay modes ip(2S) — > ~yXc2, (a) Xc2 ~^ 
2(7T + 7r _ ) and (b) Xc2 K + K~ir + Tr~ir° . The points are data; the dashed lines are the signals; the 
dotted lines are the backgrounds; and the solid lines are the sums of signal and background. 



with the exception of the 3(7r + 7r~) final state. The value listed by the PDG for B(x C 2 — > 
3(7r + 7r - )) comes from one measurement [20]. The same paper presents a measurement 
for B(xc2 2(7r + 7r~)) that leads to an expectation of 5112 ± 1334 observed events, well 
below both the PDG expectation and our measurement. The ratios of the current yields 



to those derived from Ref. |20j are 1.4 ± 0.4 and 1.8 ± 0.4 for the 2(7r + 7r _ ) and 3(7r + 7r~) 
final states, respectively. Our observed yield for the 2(7r + 7r~) final state is consistent with 
the current PDG value for B(x C 2 2(7r + 7r~)), which is determined from a 28-parameter fit 
using properties of the Xcj and ip(2S). While further measurements may clarify the B{xc2 
3(7r + 7r~)) discrepancy, we conclude that the ip(2S) — > -fXc2,Xc2 — > X study satisfactorily 
validates our yield determination procedure. 



IV. YIELD DETERMINATIONS 

Figures [3] - [5] show the measured energy distributions of the transition photon candidates 
in the rj c (2S) signal region. The photon energy resolution in this region, determined from 
the signal MC samples, is ~4.6 MeV. No significant signal is observed in any mode. The 
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TABLE II: Expected and observed yields for the tp(2S) — > 7Y C 2 study. For each mode, Bpdg is 
the value and uncertainty of B(x C 2 X) from the PDG [14|. The values in column u Npdg" 
are determined from N PDG = B(ip(2S) -> 7Xc2)<6(Xc2 -> X)N^ 2S )€, where B(ip(2S) -> 7Xc2) = 
(8.1 ± 0.4)%, N^s) = 25.9 x 10 6 , and e is the detection efficiency The values in column "N S i g " 
are fit results and the errors are statistical only. 



Channel 


Bpnr fxl(T 3 ) 


TV p nr* 


sig 


AT,,„ - Npnr 

Sly 1 U\j 




0.71 ±0.11 


262 ± 43 


294 ± 17 


32 ±46 


K+K-tt 


0.36 ± 0.09 


192 ± 49 


219 ± 17 


27 ±52 


2(vr+7r-) 


12.5 ± 1.6 


6947 ± 953 


7215 ± 119 


268 ± 960 


3(vr+7r-) 


8.7 ±1.8 


3364 ± 716 


6083 ± 113 


2719 ± 725 


K+K-iT+iT- 


10.0 ±2.6 


5226 ± 1383 


4717 ±95 


-509 ± 1386 


K+K-1T+1T-1T 






3197 ±62 




K+K-2(ir + ir-) 






2249 ± 68 










1453 ± 54 






0.56 ±0.15 


109 ± 29 


141 ± 14 


32 ±32 


K+K-j] 


< 0.4 


< 89 


51.3 ±9.1 






0.59 ±0.22 


28 ±11 


3.7 ± 5.2 


-24.3 ±12.2 



2(7r + 7r~) decay mode is the only mode in which an excess above background is present. We 
have investigated other aspects of the events in the signal region and found that this excess, 
which has a statistical significance of slightly more than 3 standard deviations, is most likely 
caused by an upward fluctuation of the background (l8j . 

For the final states that do not include an r\ decay, the signal yield upper limits are deter- 
mined by finding the value corresponding to 90% of the probability distribution determined 
from the measurement, restricted to physically allowed values. The yield measurements are 
listed in Table IHIl 

For the final states that include an t] decay, which have very low statistics, the method of 
Feldman and Cousins 21[ is used and only 90% confidence level upper limits are determined. 
The shower energy distribution is divided into two regions: a signal region (34 — 62 MeV) 
corresponding to one full width about the T] C (2S) mass [£f and a sideband region (66 — 
94 MeV). The sideband region is fitted with the shape from the background MC samples 
and the resulting normalization is used to compute the number of background events in the 
signal region, listed as Nb g in Table IIVI The number of observed events in the signal region 
is given as N b s . 

Figure |6] shows the summed shower energy distribution for the ten rj c (2S) decay channels. 
The background distribution in Fig. E] (dotted histogram) has been constructed by adding 
mode-by-mode background-only fits. It shows a visible overestimate compared to the data 
distribution in the lowest energy bins, which we attribute to the modeling of the split- 
off distribution. The alternative treatment using the two-parameter background MC fit 
described in Sec. [III] (solid histogram) reproduces the low-energy range noticeably better. 
While there may be a small excess in the signal region above the estimated background 
with the two assumed background shapes, the statistical significance is less than 3 standard 
deviations and is dependent on the background shape. 
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FIG. 3: Distributions of measured shower energy in the rj c (2S) signal region for the modes (a) 
KKir, (b) 2(7r + 7r~), and (c) K + K~ir + ir~ . The points are data; the dashed lines are the signals; 
the dotted histograms are the backgrounds; and the solid histograms are the sums of signal and 
background. 



V. SOURCES OF SYSTEMATIC UNCERTAINTY 

The sources of systematic uncertainty in our measurements of r] c (2S) branching fractions 
have been evaluated by reanalyzing the if)(2S) — > jrjci^S) and ip(2S) — > ■JXc2 candidate 
samples with alternative procedures. They are listed in Tables HTH and [TV] and described 
in detail below. All individual uncertainties not explicitly listed in the tables are combined 
in quadrature and listed as "Other." The uncertainty in the number of ip{2S) decays in 
our sample (2% [22[) also affects the branching fraction determinations. Other uncertainties 
in the detection efficiencies that have been evaluated include those associated with trigger 
decisions (1%), reconstruction of the transition photon (2%) and other particles, and PID. 

The particle reconstruction and PID uncertainties have been estimated using 281 pb _1 
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FIG. 4: Distributions of measured shower energy in the rj c (2S) signal region for the modes (a) 
3(-k + tt-), (b) K + K-tt + tt-tt°, (c) K + K~2(ir + n~), and (d) K S K ^^tt+tt' . The points are data; 
the dashed lines are the signals [not shown for the K + K~2(tt + it~) mode since its area is less than 
zero]; the dotted histograms are the backgrounds; and the solid histograms are the sums of signal 
and background. 



of data collected at the peak of the ^(3770) resonance [23J. Reconstruction uncertainties 
are determined by detecting all particles in an event except the particle being investigated, 
determining the efficiency for reconstructing the particle in data and MC simulations, and 
taking the difference as the systematic uncertainty. These studies find uncertainties of 0.3% 
per charged pion, 0.6% per charged kaon, 1.8% per K s , and 2.0% for reconstruction of the 
77 decays of 7r° and 77. The uncertainties associated with PID are determined by comparing 
the efficiency differences between data and MC simulations after applying the PID criteria. 
The uncertainties are found to be 0.25% (0.3%) per charged pion (kaon). 

The uncertainties determined from studying ip(2S) — > 7Xc2,Xc2 — > X decays described 
below are found one by one by removing the selection criterion being investigated, rede- 
termining the efficiency-corrected yield, and taking the relative difference between this and 
the nominal case as the systematic uncertainty. For J/ip suppression using the 7r + 7r~ recoil 
mass, the uncertainties for the 3(7r + 7r~), K + K~ n + n~ -k° , and K + K~2(tt + tc~) final states are 
2.4%, 2.2%, and 1.3%, respectively; they are less than 1% for other modes. For J/ijj sup- 
pression using the invariant mass of the decay products, the uncertainties for the 3(7r + 7r~) 
and K + K~2(7t + tt~) hadronic final states are 1.0% and 1.5%, respectively; they are less 
than 1% for other modes. An uncertainty of 2.4% is assigned to the requirement on the 
angle between the initial pion momentum and the candidate photon based on studies of 
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FIG. 5: Distributions of measured shower energy in the i] c (2S) signal region for the modes (a) 
7r + 7r~?7, (b) K + K~r], and (c) 7r + 7r™?/. The points are data; the solid histograms are the back- 
grounds; the solid arrows enclose the signal region; and the dashed arrows enclose the sideband 
region. 



Xc2 2(7r + 7r~) decays. An uncertainty of 0.8% is assigned to the requirement on the dis- 
tance between a transition photon candidate and the nearest track in the CC based on 
studies of Xc2 — > K + K^n + n^ decays. An uncertainty of 2.3% is assigned for suppressing 
transition photons from ip(2S) — > ^Xcj decays in selecting 7r° and i] — > 77 candidates by 
studying Xc2 — > K + K~tt + 7t~7t°. An uncertainty of 1.3% is conservatively assigned for the 
requirement on the vertex fit by taking the uncertainty from the Xc2 decay mode with the 
largest discrepancy. An uncertainty of 4.0% (2.2%) is conservatively assigned to the full 
event fit by taking the largest deviation from the Xc2 decay mode with (without) a n° or r\ 
decay. 

The uncertainties associated with the detector resolution and minimum shower energy for 
the signal region are determined by varying these parameters in the study of Xc2 decays. The 
uncertainties are assigned by recalculating the Xc2 ~^ A efficiency-corrected yield and taking 
the difference from the nominal case. The detector resolution uncertainty is determined by 
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FIG. 6: Summed distribution of measured shower energy in the rj c (2S) signal region for all i] c (2S) 
candidates decaying to light hadrons in all ten rj c (2S) decay modes. The points are data; the dotted 
histogram is the background without special treatment of split-off showers (one free parameter); 
and the solid histogram is the background separately treating the normalizations of the split-off 
and nonsplit-off showers (two free parameters). 



TABLE III: Summary of results and systematic uncertainties for rj c (2S) modes with yields de- 
termined by fitting the transition photon candidate energy distribution. The "x 2 /d.o.f." column 
lists the fit results for Figs. [3] and HI Central values and 90% confidence level upper limits are 
provided for the product branching fraction B\B 2 = i3(ip(2S) — > "yr] c (2S))B(r] c (2S) — > X). The 
columns under "Systematic uncertainties" correspond to (A) AM criterion, (B) background pa- 
rameterization, (C) M[r/ c (2S)] uncertainty, (D) signal region, and (E) nonresonant background 
component, as described in the text. Only statistical uncertainties are included in the N S i g results, 
while statistical and systematic uncertainties are included in the B\B 2 results. We have assumed 
r[7fe(2S)] = 14 MeV/c 2 " 



Channel 



N 



.sig 



X 2 /d.o.f. Systematic uncertainties (%) 

A B C D E Other Total 



B X B 2 (xKT 6 ) 

(90% C.L. 



KKn 




11.75$ 


9.7/14 


15.2 


16.9 


12.2 


13.3 


8.2 


5.6 


30.7 


5.9^5 ± 1-8 


<14.5 


2(7T+7T 


") 


47.9±if;B 


14.1/14 


2.5 


6.1 


11.8 


11.7 


4.4 


5.0 


19.1 


9.0±|;| ± 1.7 


<14.6 


3(7T + 7T 


") 


10 1+ 18 - 1 

lu - i -17.6 


11.2/14 


16.6 


20.4 


14.6 


12.1 


3.6 


5.1 


33.0 


2.7±t| ± 0.9 


<13.2 


K+K~ 




19 c+15.8 
-^•0-15.6 


9.2/14 


7.7 


32.6 


7.1 


13.8 


4.4 


4.5 


37.2 


2.5±|;1±0.9 


<9.6 


K+K- 


7T + 7r~7r° 


07 r+21.3 

°' -°-20.8 


13.8/14 


22.7 


15.2 


29.7 


24.5 


0.9 


7.2 


47.8 


16.7±H ± 8.0 


<43.0 


K+K- 


2(vr+7r-) 


n q+12.6 

u -°-12.2 


13.2/14 


0.8 


6.8 


14.8 


11.0 


4.3 


5.3 


20.8 


-0.lt|; 7 ±0.1 


<9.7 


K°K± 




12.9+f;f 


11.9/14 


13.2 


17.8 


16.3 


5.6 


5.2 


5.0 


29.0 


6.4+|; 7 ± 1.8 


<15.2 
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TABLE IV: Summary of results and systematic uncertainties for rj c (2S) modes with yields deter- 
mined by sideband subtraction. The product branching fraction is defined as B1B2 = B(ip(2S) — > 
'yr] c (2S))B(ri c (2S) — > X). The columns under "Systematic uncertainties" correspond to (A) AM 
criterion, (B) background parameterization, and (C) M[rj c (2S)} uncertainty, as described in the 
text. Statistical and systematic uncertainties are included in the -B1-B2 results. We have assumed 
r[?y c (2S)] = 14 MeV/c 2 Q. 





N obs 1 N bg 


N sig /e 


Systematic uncertainties (%) 


B^B 2 (xlO" 6 ) 


Channel 




(90% C.L.) 


ABC Other Total 


(90% C.L.) 




4 / 4.3 


<75.4 


6.1 46.5 3.0 8.1 47.7 


<4.3 




8 / 6.5 


<115.7 


9.2 29.8 3.9 7.6 32.3 


<5.9 




2 / 1.8 


<287.9 


8.9 24.5 3.7 7.4 27.3 


<14.2 



individually varying the width of the core Gaussian and the transition point between the 
core Gaussian and the power law tail of the Crystal Ball function by 1 standard deviation, 
resulting in a total uncertainty of 1.0%. An uncertainty of 3.2% is assigned to the minimum 
photon energy requirement by increasing the lower bound from 90 to 110 MeV, corresponding 
to the same difference (18 MeV) between the minimum energy and the r] c (2S) mass in the 
T] C (2S) signal region. 

The uncertainties associated with selecting 77 — > 7r + 7r~7r° and 7/ — > 7r + 7r~77 decays based 
on the invariant mass of the decay products and the 77 recoil mass used for i()(2S) — > rjj/ip 
suppression are determined by studying ip(2S) —> r yXc2,Xc2 — > X decays. The individual 
uncertainties are assessed by varying the respective mass range to double the detection inef- 
ficiency, redetermining the Xc2 — > X efficiency-corrected yield, and assigning the difference 
between this and the nominal result as the systematic uncertainty. An uncertainty of 1.3% 
is assigned for the 7r + 7r~7r° invariant mass selection range in the 7r + 7r~77 and K + K~r] decay 
modes. An uncertainty of 1.2% is assigned to the 7r + 7r~77 invariant mass selection range for 
the 7r + 7r~?/ decay mode. Uncertainties of 3.7% and 0.2% are assigned to the 77 recoil mass 
suppression range for the 77 — > 77 and 77 — > 7r + 7r~7r° decays, respectively, comprising the 
7r + 7r~77 decay mode. 

The largest systematic uncertainties in all r] c (2S) decay modes arise from the AM se- 
lection criterion, the parameterization of the background shape, and the uncertainty in the 
i] c (2S) mass, which are listed as separate entries in Tables [TTT1 and HVl They are estimated by 
varying the criterion being investigated, redetermining the efficiency-corrected yield upper 
limit in the rj c (2S) signal region, and assigning the difference between this and the nominal 
result as the systematic uncertainty. The uncertainty associated with the AM criterion is 
determined by removing the cut. The uncertainty arising from the background modeling 
for the decay modes without an r? decay is determined by replacing the background deter- 
mined from the background MC samples with a first-order polynomial. The uncertainty 
arising from the background modeling for the decay modes with an rj decay is determined 
by lowering the overall background yield in the signal region by 1 standard deviation based 
on the data yield in the sideband region. The effect of the uncertainty of the r] c (2S) mass is 
determined by 1 standard deviation variations of the mass, M[?7 C (2S')] = 3638±4 MeV/c 2 , in 
the fits of the measured shower energy distributions and the determination of the detection 
efficiencies, with the larger discrepancy from the two cases being assigned as the systematic 
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uncertainty. 

Additional systematic uncertainties in rj c (2S) decay modes without an 77 decay arise 
from the signal region range and the nonresonant component of the background, which 
are listed as separate entries in Table IHIl The uncertainty associated with the maximum 
boundary of the signal region is assessed by varying the boundary by 8 MeV, with the larger 
deviation from nominal of the two cases being assigned as the systematic uncertainty. The 
uncertainty associated with the minimum boundary of the signal region described above is 
combined in quadrature to obtain the signal region uncertainties listed in Table IHIl The 
uncertainty arising from the nonresonant component of the background was investigated by 
determining the ratio of event yields in the 20.6 pb" 1 of off-resonance data collected at \/s 
= 3.67 GeV with the yields from a 5 times luminosity continuum MC sample generated at 
the same center-of-mass energy. The same event selection criteria were applied as for the 
r] c (2S) signal search with the exception that AM was redefined as AM = 3.67 GeV — M inv . 
The uncertainty was assessed by repeating the fits of the r] c (2S) signal region but with the 
continuum MC component of the background fixed to the ratio found in the off-resonance 
data sample study. The difference between the rj c (2S) signal yields from this and the nominal 
result was assigned as the systematic uncertainty. No systematic uncertainties were applied 
to the 77 decay modes due to the small amount of nonresonant background. 

Tables IHIl and IIVI summarize the total systematic uncertainties. The individual uncer- 
tainties are treated as uncorrelated and are combined in quadrature to obtain the overall 
systematic uncertainties in the product branching fraction upper limits. The total uncer- 
tainty for composite decay modes is determined by weighting the total systematic uncertainty 
of each constituent decay mode by its branching fraction. 

In addition to these sources of error, the partial width for a direct Ml radiative transition 
between ip(2S) and rj c (2S) is related to the matrix element governing the spin-flip transition 
I and the energy of the transition photon E 1 by 

T^(2S)^ iVc(2S)] oc Ep 2 . (3) 

This implies a signal shape given by a Breit-Wigner times function, rather than the Breit- 
Wigner function that was used for our fits. We studied the effect of using this modified 
signal shape on the product branching fractions with the following procedure. For the 
determination of yields in i] c (2S) modes without an 77 decay, the signal regions were fitted 
with Breit-Wigner times E^ functions convoluted with Crystal Ball detector resolution signal 
shapes and the histogram backgrounds used in the nominal results. The nominal yields were 
used for modes with an rj decay. The efficiency for each mode was determined by applying 
the nominal event selection criteria to signal MC samples generated with signal shapes 
that were Breit-Wigner distributions multiplied by E^. Since there is no obvious choice 
of damping function as there was for the ground state resonance jsj, we use an arbitrary 
cutoff on the maximum allowed photon energy. With this procedure, we find deviations in 
the product branching fractions that are on the order of, and in some cases greater than, 
the other uncertainties. Because of the arbitrary cutoff, it is difficult to assign a systematic 
uncertainty to this effect, and we have chosen not to include it. 

VI. BRANCHING FRACTION RESULTS 

The upper limits on the number of signal events are used to set upper limits on the 
product branching fractions. The product branching fraction for each rj c (2S) decay mode is 
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determined by Eq. fl2]), where N sig is the number of signal events and N^s) — 25.9 x 10 6 . 
Because no statistically significant signals are observed in any of our ten decay channels, 
we use the efficiency-corrected yields to set upper limits on the product branching fractions. 
Systematic uncertainties are determined and combined as described in Sec. [V] and added 
to the statistically calculated product branching fraction upper limits. The final results are 
given in Tables HTT1 and HVl 

Upper limits for the product branching fraction are also determined as a function of 
the T] C (2S) full width, for which the current world average is T[q c {2S)} = 14 ± 7 MeV/c 2 
jsj. Separate signal MC samples with r[?7 c (2S')] = 7 and 21 MeV/c 2 were generated in the 
same manner as the nominal MC samples. The measured shower energy distributions are 
fitted in the same manner as for the standard yield determination procedure, but with the 
resolution functions determined from these MC samples and the full width of the signal 
shape adjusted to match the full width being investigated. The linear extrapolation of the 
product branching fraction as a function of r[r/ c (2S')] is listed in Table |V] for each r] c (2S) 
decay mode. 

TABLE V: Summary of product branching fraction results as a function of T[r] c (2S)}. The y- 
intercept and slope parameters a and b are defined by B{ip{2S) — > jr] c (2S))B(r] c (2S) — > X) < 
a + b x r[?7 c (2S')]. Statistical and systematic uncertainties are included in these results. 





a 


b 


Channel 


(io- 6 ) 


(10" 6 c 2 /MeV) 


KKtt 


6.6 


0.56 


2(7T+7T-) 


6.5 


0.58 


3(7r+vr-) 


4.0 


0.74 


K+K--K+H- 


3.1 


0.50 


K+K-ir+Tr-ir 


15.5 


2.09 


K + K-2(ir + ir-) 


5.2 


0.34 




8.5 


0.51 




2.9 


0.09 


K+K-j] 


4.1 


0.13 




10.0 


0.31 



VII. SEARCH FOR t) c (2S) 7t+7t-t] c (1S) 

In addition to searching for t] c (2S) decays to light hadrons, a search for the decay pro- 
cess ip{2S) — > r yr] c (2S),r] c (2S) — > n + n^ri c (lS) is also performed. The four hadronic final 
states of 7r + 7r~ ± 7T T: ) , tt+tt- (K + T^tt ), 7r+7r-[2(7r+7r-)], and tt + 7t-(K + K-tt + 7t-) plus 
a candidate transition photon are used for this study. The selection criteria are the same as 
described above, except that the J/ip rejection criterion based on the ir + n~ recoil mass is 
removed and we require the hadronic decay products not associated with the dipion transi- 
tion to be within 40 MeV/c 2 of M[rj c (lS)]. Information from the n + ir~ recoil mass is not 
used since the distribution is broadened by the intrinsic widths of the ^(l^) and rj c (2S). 
Figure [7|(a) shows the invariant mass of the rj c (lS) candidates. 
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FIG. 7: (a) Hadronic invariant mass of rj c (lS) candidates. The points are data; and the solid 
histogram is signal MC, arbitrarily normalized for clarity. Selected rj c (lS) candidates are enclosed 
by the arrows. All other event selection criteria have been applied, (b) Measured shower energy 
distribution for candidates in the rj c (2S) signal region after applying the selection criterion for the 
ift(2S) — > r yrj c (2S),rj c (2S) — > TT + TT~r] c (lS) decay. The points are data; the dotted histogram is the 
background; and the solid histogram is the sum of signal and background. The signal is not shown 
since its area is less than zero. 

Figure [T^b) shows the measured shower energy distribution after applying the T] C (2S) — > 
7T + 7T~r] c (lS) selection criterion. No evidence of a signal is observed. The xVd.o.f. of the 
fit is 19.5/14. The number of signal events is N sig = — 5.lig°g 7 , corresponding to an upper 
limit of N sig < 14.8 (90% C.L.). 

The upper limit on the product branching fraction is determined by 

N ■ 

Nf(2S)[Ei€i x Bi{r] c (lS))\ 

where, for a given final state i, is the detection efficiency and Bi(rj c (lS)) is the branching 
fraction for the r) c (lS) decay. The signal efficiency for each r) c (lS) decay mode is determined 
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from signal MC samples. The decay ip(2S) — > 7?7 C (2S') is generated in the same manner as 
described in Sec. HO The decay rj c (2S) — > ir + 7r~r] c (lS) is generated according to phase 
space with the r] c {lS) parameters M[r] c (lS)] = 2979.8 ± 1.2 MeV/c 2 and V[r] c (lS)} = 26.5 ± 
3.5MeV/c 2 [14]. The individual values used to determine £i e; xBi(r] c ( IS)) = (0.50±0.07)% 
are listed in Table IVT1 



14 MeV/c 2 [8j, ll4j. The column 



TABLE VI: Efficiency and submode branching fraction information for the ip(2S) 
'yr) c (2S),r) c (2S) — > tt + tt^i] c (1S) study. We assume T[r] c (2S) 
B(r] c (lS)) lists the branching fractions for the rj c (lS) decay 
submode branching fractions. 



14j | . The listed efficiencies include 



r) c (lS) Decay Mode 




B( Vc (lS)) (%) 


e x B( Vc (lS)) (%) 




6.64 ±0.14 


2.3 ±0.4 


0.148 ±0.026 


K+K-1T 


7.33 ±0.18 


1.17 ±0.20 


0.086 ±0.015 


2(7r+vr-) 


10.99 ±0.21 


1.2 ±0.3 


0.13 ±0.03 


K+K-iT+iT- 


8.88 ±0.20 


1.5 ±0.6 


0.13 ±0.05 



Several sources of systematic uncertainty in the r) c (2S) — > ir + ir~r) c (lS) measurement 
have been evaluated. The uncertainty due to the invariant mass range used to select r) c (lS) 
candidates is determined by tightening the mass range to double the detection inefficiency, 
redetermining the product branching fraction upper limit, and assigning the relative dif- 
ference between this and the nominal result as the systematic uncertainty. The effect 
of the uncertainty of the 77 C (1S) full width is determined by generating separate signal 
MC samples with rf^lS 1 )] = 23 and 30 MeV, i.e., 1 standard deviation variations of 
r[?7 c (lS')] = 26.5 ± 3.5 MeV/c 2 [l4|, to redetermine the detection efficiencies and repeating 
the yield determination procedure with the resolution functions determined from these MC 
samples. The uncertainties associated with the AM criterion, background parameterization, 
M[i] c (2S)} uncertainty, signal region, and nonresonant background component are evaluated 
by performing the procedures described in Sec. [V] The systematic uncertainties associ- 
ated with the number of ip(2S) decays, trigger efficiency, particle reconstruction, PID, full 
event and vertex fitting, and the suppression of transitions to other charmonium states for 
these specific final states are assigned as described in Sec. |V] Table [VTT1 lists the individual 
contributions to the total systematic uncertainty of the product branching fraction. 

The product branching fraction for ^(2,5) — > r yrj c {2S) ) r q c {2S) — > 7i + 7i^ri c (lS), assuming 
r[r] c (25)] = 14 MeV/c 2 and including statistical and systematic uncertainties, is B(ip(2S) — > 
jr] c (2S))B(r] c (2S) n+n-^lS)) = (-0.39±° o f 6 ± 0.18) x 10" 4 < 1.7 x 10~ 4 (90% C.L.). 
Expressed as a function of T[r] c (2S)], the upper limit is {(48) + (9.2 c 2 /MeV) x T[r] c (2S)]} x 
10~ 6 . 



VIII. SUMMARY AND CONCLUSIONS 



In summary, we do not observe the transition ip(2S) — > 7?7 C (2S') with any of the ten 
exclusive r] c (2S) decays to light hadrons. We also do not observe evidence for the decay 
process ip{2S) — > r )Tj c {2S) ) rj c {2S) — > it + it~ri c {lS). Our original objectives for measuring the 
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TABLE VII: Sources of systematic uncertainties for the product branching fraction B(ip(2S) —> 
7 t] c (2S))B(t] c (2S) — > ir + ir~7] c (lS)). The total systematic uncertainty is determined by combining 
the individual contributions in quadrature. 



Source 


Uncertainty (%) 


Nonresonant background 


28.6 


M[ Vc (2S)} 


22.0 


Signal region 


20.7 


B(Vc(lS)) 


14.0 


M[ Vc (lS)] 


11.1 


AM 


7.9 


nvc(is)} 


4.1 


Background parameterization 


1.5 


Other 


5.5 


Total 


46.5 



properties of the r] c (2S) cannot be achieved with this data sample, and only upper limits 
for the product branching fractions are obtained. 

The BaBar Collaboration recently reported a branching fraction of B(r] c (2S) — > KKtt) = 
(1.9±0.4(stat) ±l.l(syst))% where the systematic uncertainty is dominated by the in- 
clusive measurement of B ± — > K ± r] c (2S) [25}. Using the central value of B(t] c (2S) — >■ KKtt) 
and our 90% confidence level upper limit of B(ip(2S) -»■ -fr] c {2S))B{rj c (2S) ->• KKtt) < 
14.5 x 10 -6 leads to B(ip(2S) — > , yr] c (2S)) < 7.6 x 10~ 4 , which is larger than the phenomeno- 
logical prediction of B{tp{2S) -»■ 7^(25)) = (3.9 ± 1.1) x 10" 4 . 

The B(rj c (2S) — > KKtt) measurement can also be used to determine upper limits of t] c (2S) 
hadronic decays based on published rj c (2S) searches. The two-photon fusion result reported 
by the CLEO Collaboration ^ and the B(t] c (2S) — > KKtt) measurement lead to a two- 
photon partial width of r 77 [?7 c (2S')] = 4.8 ±3.7 keV. Using this value of r 77 [?7 c (2S')] with the 
recent two-photon fusion upper limits for r] c (2S) production from the Belle Collaboration 
0, we find B(r) c (2S) -> 2(7r+7r-)) < 0.14% and £(^(25) -> K+K'tt+tt') < 0.10% 
(90% C.L.). These upper limits are an order of magnitude smaller than the branching 
fractions obtained by assuming that the partial widths for T] C (2S) decays are the same as for 
Tj c {13), i.e., B(r) c {2S) -> 2(7r+7r-)) = (2.3 ± 0.6 ± 1.2)% and B{r] c {2S) ->■ K+K-tt+tt-) = 
(2.9 ±1.1 ±1.5)%, where the first error is the uncertainty from the rj c (lS) branching fraction 
and the second error is the uncertainty from r[?y c (2S')] (8). 
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